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We examine the fate of supersymmetric fiat directions. We argue that the non-perturbative decay 
of the fiat direction via preheating is an unlikely event. In order to address this issue, first we 
identify the physical degrees of freedom and their masses in presence of a large flat direction VEV 
(Vacuum Expectation Value). We explicitly show that the (complex) flat direction and its fermionic 
partner are the only light physical fields in the spectrum. If the flat direction VEV is much larger 
than the weak scale, and it has a rotational motion, there will be no resonant particle production at 
all. The case of multiple fiat directions is more involved. We illustrate that in many cases of physical 
interest, the situation becomes effectively the same as that of a single flat direction, or collection of 
independent single directions. In such cases preheating is not relevant. In an absence of a fast non- 
perturbative decay, the flat direction survives long enough to affect thermalization in supersymmetric 
models as described in hep-ph/0505050 and hep-ph/05I2227. It can also "terminate" an early stage 



of non-perturbative inflaton decay as discussed in hep-ph/0603244 



The Minimal Supersymmetric Standard Model 
(MSSM) has a large number of flat directions, along 
which the scalar potential vanishes. The flat directions 
are parameterized by gauge-invariant monomials [l|, 0] , 
for a review see [§]■ Recently it has been pointed out 
that these flat directions can play an important role in 
thermal history of the Universe within supersymmetry 
(SUSY) [J, Q . We remind that the thermalization rate 
depends on 2 <-» 2 and 2^3 scatterings of the inflaton 
decay products. The former lead to kinetic equilibrium 
but do not change the total number of particles. How- 
ever, after the last stage of (perturbative) inflaton decay, 
the reheat plasma is far from chemical equilibrium as 
the number of particles is much less than that in full 
thermal equilibrium [4]. As pointed out in [6], 2 — > 3 
scatterings with gauge boson exchange in the i-channel 
play the main role in establishing chemical equilibrium. 
For massless gauge bosons, they are infrared divergent. 
In a non-supersymmetric set up the plasma mass of 
gauge bosons provides a natural infrared cut-off which 
regulates these processes(for details see Q). Since right 
after the inflaton decay the reheat plasma is dilute, the 
resulting plasma mass is small. Thus 2^3 scatterings 
are enhanced and additional particles are efficiently 
produced. In consequence, full thermal equilibrium is 
established quickly in this case [fj. 

However the situation will be completely different 
within supersymmetry because of the flat directions. 
During inflation these flat directions develop large Vac- 
uum Expectation Value (VEV), see for details [3]. A 
large VEV induces large SUSY preserving masses to 
gauge bosons and gauginos via the Higgs mechanism [4| 1 . 



1 There are ways to suppress the flat direction VEV, if the MSSM 
flat direction is heavy during inflation then the flat direction set- 
tles down to the minimum of its potential. However this happens 
for a very special case of the Kahler potential, see for details 



The 2 — ► 3 scatterings mediated by the gauge bosons 
are therefore strongly suppressed and, consequently, ther- 
malization will be very slow [J, Q . We also showed that 
for a Standard Model (SM) singlet inflaton the inflaton 
decay is typically perturbative within SUSY 4, 8] 2 . 

Delayed thermalization has important outcomes: the 
Universe undergoes a quasi-thermal phase for a long pe- 
riod after the inflaton has completely decayed, during 
which the reheat plasma is dilute and far from chem- 
ical equilibrium [j, After the flat direction starts 
oscillating, its VEV and its induced mass for gauge 
bosons/gauginos is redshifted by the Hubble expansion. 
Gradually thermalization rate increases and full thermal 
equilibrium is established shortly after this rate becomes 
comparable to the expansion rate of the Universe. A 
full thermalization typically happens much later than the 
inflaton decay, leading to a reheat temperature (much) 
lower than naively expected Q. The exact value of the 
reheat temperature depends on the nature of the flat di- 
rection and its initial VEV and, typically it is in the range 
T R - 10 4 - 10 7 GeV [|. 

The notable off-shoots of a delayed thermalization and 
a low reheat temperature are 3 : 

• Thermal gravitino production is suppressed. SUSY 
provides a natural solution to the infamous thermal 
gravitino problem 0, [5[ . Note that the abundance 
of non-thermal gravitinos is not affected [14 . fig . 



2 In almost all models of inflation, the inflaton is treated as a 
SM gauge singlet. Very recently we constructed models of in- 
flation which are gauge invariant combinations of squarks and 
sleptons [TJ|, and with right-handed sneutrinos |10| . The model 
predictions are trustable and robust (see for other gauge invari- 
ant inflatons We emphasize that this novelty is solely due 
to the fact that the inflaton is charged under the MSSM gauge 
group. 

3 A late thermalization holds the key to a successful MSSM cur- 
vaton mechanism [l2j|. see also [l3| . 



2 



but they do not give rise to a large gravitino abun- 
dance any way. 

• Thermal leptogenesis within MSSM does not work 
in most of the parameter space [3, Q • 

The above conclusions hold if the flat direction VEV 
remains large for a sufficiently long time after the onset 
of its oscillations. This is the case if the flat direction 
decays perturbatively. Note that the flat direction VEV 
induces a large mass to all fields which are coupled to 
it. Perturbative one-particle decay to these fields is only 
possible when their mass, hence the flat direction VEV, 
is less than the flat direction mass. If the initial VEV 
is much larger than the weak scale, this happens much 
later than the oscillations start. 

Recently it has been pointed out in [l6[ that, the flat 
direction oscillations can decay very rapidly (i.e., within 
a few oscillations) via preheating [17|, LL8[ ■ This will have 
important implications for cosmology. For example the 
Universe will quickly thermalize even in the presence of 
a flat direction with a large VEV. Here we examine this 
possibility in some detail and argue that the decay of 
MSSM flat directions via preheating is unlikely in cases 
of physical interest. 

First of all, we identify the physical degrees of freedom 
and the flat direction masses when a single flat direction 
has a non-zero VEV. We show that the flat direction and 
its fermionic partner are the only light fields in the rel- 
evant spectrum. As a result, the flat direction motion, 
so long as it has a rotational component, results in time 
variation of the mass eigenstates and eigenvalues of su- 
perheavy fields. This, however, is not significant at all and 
does not lead to resonant particle production [2(| Hl| . 

Next, we consider the case with more than one flat 
directions. Two or more flat directions which are simul- 
taneously D- and F-flat can develop large VEVs at the 
same time. Within MSSM the only directions which are 
completely independent are the two represented by the 
LLe and udd monomials. In this case both the flat direc- 
tions evolve independently and their dynamics is equiva- 
lent to a single flat direction, and hence there will be no 
preheating. We also emphasize that preheating will not 
be relevant when two flat directions have very different 
VEVs. This happens in many cases since the directions 
under consideration are lifted by different F-terms which 
become important at different field values. 

The second possibility is multiple flat directions associ- 
ated with a gauge- invariant polynomial (when all families 
of leptons and quark are taken into account). This is a 
new case and very involved, even for a simple example 
of H.„L monomial, whose dynamics was first described in 
Ref. [22|, but preheating effects were not discussed. In 
a general case, however, the motion of the flat direction 
VEV is not known yet. One cannot embark on the is- 
sue of particle creation before knowing the details of the 
flat direction motion for such a scenario. We will briefly 
sketch this scenario and we think preheating is unimpor- 
tant for practical purposes. Moreover, due to F-term 



lifting, in many cases the flat manifold has complex di- 
mension one at very large field values. In such cases the 
situation becomes similar to that of a single flat direction. 

The central message is: preheating of an MSSM fiat 
direction is unimportant for physically interesting cases, 
thus, the flat direction lives long enough to affect reheat- 
ing and thermalization. 



I. SPONTANEOUS SYMMETRY BREAKING 
AND PHYSICAL DEGREES OF FREEDOM 

A crucial point is to identify the physical degrees of 
freedom and their mass spectrum in presence of a non- 
zero flat direction VEV 4 . 

Let us consider the simplest flat direction, which in- 
cludes only two fields: H u Hd- This is also familiar from 
the electroweak symmetry breaking in MSSM. A clear 
and detailed discussion is given in [23|]. Here we briefly 
recount the situation. 

One can always rotate the field configuration to a ba- 
sis where, up to an overall phase, If„ = H\ — and 
H\ = Hi = i^o/v^-Here superscripts denote the weak 
isospin components of the Higgs doublets. In this basis 
the complex scalar field is defined by: 



V2 



(1) 



represents a flat direction. Its VEV breaks the SU(2)w X 
U(1)y down to U(l) em ( m exactly the same fashion as in 
the electroweak vacuum) . The W and Z gauge bosons 
then obtain masses my/ , mz ~ g4>o from their couplings 
to the Higgs fields via covariant derivatives (g denotes a 
general gauge coupling). 
We also have: 



Xi 



{Hi -H\) 
V2 



and 



X2 



(H u 



Hi) 



V2 



(2) 



Then \2 and Xi R (-^ an d I denote the real and imagi- 
nary parts of a complex scalar field respectively) acquire 
masses equal to mw and mz, respectively, through the 
D— term part of the scalar potential. Note that 



X3 



(Hj -Hi) 
V2 



(3) 



and Xi, i are the three Goldstone bosons, which are eaten 
by the gauge fields via the Higgs mechanism. Therefore, 
out of 8 real degrees of freedom in the two Higgs doublets, 
there are only two light physical fields: cpn, tpj. They are 



4 The first version of Ref. Il6t did not count correctly the light de- 
grees of freedom. Then it was claimed that even a single MSSM 
flat direction can decay non-perturbatively. However, as we ar- 
gue below, this is not the case. 
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exactly massless when SUSY is not broken (and there is 
no /x term either). 

An important point is that the masses induced by the 
flat direction VEV are SUSY conserving. One therefore 
finds the same mass spectrum in the fermionic sector. 
More specifically, the Higgsino fields H\ and f/J are 
paired with the Winos, while — H^)/y/2 is paired 
with the Zino to acquire masses equal to mw and mz, 
respectively, through the gaugino-gauge-Higgsino inter- 
action terms. The fermionic partner of the flat direction 
(11% + HY)/y/2 remains massless (note that the photon 
and photino are also massless, but not relevant for our 
discussion). 

In reality, supersymmetry is broken and, ip obtains 
a mass m v ~ O(TeV) from soft SUSY breaking term 
(the same is true for the gauginos). However, for g(po ^> 
O(TeV), which is the situation relevant to the early Uni- 
verse, the mass spectrum is hierarchical: Xi, R, X2i and 
gauge fields (plus their fermionic partners) are super- 
heavy. 

In a general case the total number of light scalars, 
Nught , is given by: 



Nlight = Ntotal — (2 X Nbroken), 



(4) 



where N total is the total number of scalar degrees of free- 
dom, and N oro k en is the number of spontaneously broken 
symmetries. Note that the factor 2 counts for the num- 
ber of eaten Goldstone bosons plus the number of de- 
grees of freedom which have obtained large masses equal 
to those of the gauge bosons. In the case of H u Hd direc- 
tion, Eq. (J) reads: 2 = 8 - (2 x 3). 



II. PARTICLE PRODUCTION FROM A SINGLE 
FLAT DIRECTION 

A generic MSSM flat direction has .F-term and D- 
term couplings to other fields. The dynamics and non- 
perturbative particle creation for a flat direction with F- 
term couplings are unimportant. This has been clearly 
illustrated in Refs. [U [2J| , here we will not discuss this 
case. We will only concentrate on the £>-terms. 

For the H u Hd flat direction, supersymmetric D-term 
contribution lead to: 

Vd = (ai5i + "2^2) [<PRXi, R + P1X1, if + 
a-idl [<PRX2, R + V1X2, if + a 2 g\ [<finX3, R + ¥1X3, if , 

(5) 

with ax, &i ~ ©(lCU 1 ). Here g\ and g 2 are gauge cou- 
plings of U(l)y and SU(2)w For a general monomial 
the D-term contribution will be similar when the scalars 
are decomposed to the flat direction tp and orthogonal 
fields x- 

Let us now consider the case where the flat direction 
VEV is oscillating. The oscillations start when the Hub- 
ble expansion rate drops below the flat direction mass, 



i.e., when H(t) ~ m v . Note that the motion is due to soft 
supersymmetry breaking terms and the _F-terms as 
Vd identically vanishes for a flat direction. The soft mass 
and superpotential contribution depends on the modulus 
of flat direction \ip\. On the other hand, the "A"-terms 
result in a phase-dependent contribution, which exerts 
a torque and generates a rotation in the if plane. As 
a generic feature the flat direction motion has a major 
rotational component 



A. Rotation on a circle 

If the real and imaginary parts of the flat direction 
oscillate with the same amplitude, tpo, but with a tt/2 
phase difference, the trajectory will be a circle in the 
ip plane: ip — (poexp(im v t) . This is not what exactly 
happens in the early Universe, but as we just pointed out, 
this is a very good approximation of a generic situation. 

The first term on the right-hand side of Eq. (O, which 
arises from the diagonal part of SU (2)w and U(1)y, gives 
rise to 5 : 

V ~ (gf+gt) pi [cos(m v t)xi, r + sin(m v £)xi, if ■ (6) 

At a first glance it seems that there are two scalar mass 
eigenstates: 



ui(t) = cos(m v t)xi, r + sin(m v ,t)xi, / , 
u> 2 (t) = cos(m v i)xi, / - sm(m 9 t)xi, r ■ 



(7) 



where oj\ has a mass ~ gfoi hence superheavy, and u> 2 
is light. Note that the mass eigenvalues are constant in 
time (despite rotation of tp) but the mass eigenstates are 
evolving. 

In general one can picture an instantaneous mass eigen- 
state by a vector, to(t), whose magnitude represents the 
corresponding mass eigenvalue. In the case at hand the 
mass eigenstates are two-dimensional vectors in the xi 
plane. One can then expect non-perturbative production 
of particles whenever the adiabaticity is violated in the 
evolution of to, i.e., such that 



du(t) 



(II 



(8) 



It can be shown from Eq. ([7]) that this is the case 
at all moments of time for the light mass eigenstate 
Q 2 {t). One might then conclude that quanta of the light 
mode are copiously produced via parametric resonance, 
thus, leading to a quick decay of the flat direction 
oscillations [16|. 



5 Contributions from non-diagonal parts of SU(2)yy lead to a sim- 
ilar contribution with xi being replaced by X2, X3- 
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However, a closer inspection shows that to^it) actually 
is one of the three Goldstone bosons in the background 
of rotating ip, see Eq. and the subsequent discussions. 

Hence UJ2 (t) is not a physical degree of freedom as it is 
eaten by the Z gauge boson. Indeed, as we elaborated, 
there are no physical light fields except the rotating flat 
direction. Obviously the mass eigenstates and eigenval- 
ues of the flat direction are not affected by its rotation 
and therefore they are constant in time, see Eq. (fTJ). 

Further note that the rotation only results in a time 
variation in the mass eigenstates of the heavy fields, 
and the corresponding eigenvalues (which are constant 
in time) are much larger than the frequency of rotation. 

At the onset of the flat direction oscillations, we typ- 
ically have, gipo > 10 s m ip . This implies that the adia- 
baticity condition always holds, see Eq. ([5]), and hence 
there will be no resonant particle production 6 . 

In an absence of particle production, the flat direction 
VEV is only subject to redshifting by the Hubble expan- 
sion. In a Universe dominated by relativistic particles, 
not necessarily thermalized, we have (ip) oc H 3 / 4 . There- 
fore, the heavy fields become continuously lighter. 

The resonant particle production, due to time- 
variation of the mass eigenstates, if at all becomes impor- 
tant, will only be relevant when g(ip) ~ m v . However, 
this will happen much later than the onset of the flat 
direction oscillations, indeed, we find: 



H, 



preheat 



( m v \ 
V.9<A)/ 



4/3 



m v < 10 



m 



(9) 



where 



O(TeV). This implies that the flat di- 



rection VEV survives a very long time, which is suf- 
ficient to affect thermalization along the lines pursued 
in Refs. [3, [f|. In particular, Eq. @) can be trans- 
lated into an upper bound on the reheat temperature: 
T R ~ (HpreheatMp) 1 / 2 < 10 6 GeV. This is compatible 
even with the most stringent limit on the reheat temper- 
ature, which is obtained for unstable gravitinos with a 
dominant hadronic decay mode (l9j . 



B. A General Rotation 

For a general rotation we will have: tpn = tpo cos(m v t) 
and ipi = aipo cos{m v t + 9), where a < 1 is a positive 
number and 9 is an 0(1) phase. 

In this case, both the mass eigenstates and mass eigen- 
values of the heavy fields change in time. The mass eigen- 
values, given by ~ g\f\, he in the range ~ [agip^, g(fo\- 



Particle production will be insignificant so long as the 
time evolution of the system is adiabatic. The adiabatic- 
ity condition for mass eigenvalues is violated if: 



d \f\ ^ 2| |2 



(10) 



This happens when g\(p\ < a(gipom v ) 1 / 2 , and provided 

that a < {m^/ gipo) 1 / 2 0, HH- For a typical initial con- 
dition, where gipo > lO 8 ™^, this requires that, a < 10~ 4 . 
As mentioned earlier, in a generic situation we have 
a -0(1) 0. 

For such small values of "a" , the flat direction motion 
will be effectively one-dimensional: ipn ipocos{m v t) 
and ifi « 0. In this case, _D-terms result in an interaction 
term g 2 (p 2 R XR.> see Eq. ©, which leads to a fast decay 
of flat direction oscillations via the standard picture of 
preheating However, the challenge is to seek, a < 

10~ 4 , in a physical models. 

We conclude that in general; the flat direction oscil- 
lations survive a very long time. A fast decay can only 
happen in the exceptional case when the flat direction mo- 
tion is effectively one- dimensional. 



III. TWO OR MORE FLAT DIRECTIONS 

A necessary (but not sufficient) condition for a quick 
decay of an MSSM flat direction via preheating is that; 
there must exist light scalar degrees of freedom whose 
mass eigenstates and/or eigenvalues are evolving in time 
in a non-adiabatic fashion. As we showed this does not 
happen for a single flat direction. Here we examine the 
situation where there are multiple flat directions. Wc 
discuss two possible cases separately. 



A. Two directions represented by different 
monomials 

Two directions, which are simultaneously D- and F- 
flat, can independently acquire large VEVs during infla- 
tion. An an important example, there exists two flat di- 
rections represented by the udd and LLe monomials, re- 
spectively. Note that u and d have no Yukawa couplings 
to L and e. This guarantees that they are simultaneously 
_F-flat. The D— flatness imposes that the two d's in udd, 
and the two Us in LLe are from different families. Also 
i^-flatness requires that e belong to a third lepton fam- 
ily. Without any loss of generality, let us consider two 
flat directions, parameterized as: 



n/3' 



(up to a phase) . 



Note that for a circular rotation, the Lagrangian for the heavy 
fields is the same as that in a static case up to corrections of 
order (m v / gtpo), coming form time-derivatives of tp. Therefore, 
one can reliably use the same mass eigenstates and eigenvalues 
as that of the static case. 



L\ = L\ = e3 = (up to a phase). 



(11) 



Here the upper indices indicate color (for u and d) and 
weak isospin (for L). The lower indices denote the family. 
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Note that u can be from any of the quark families, and 
we have chosen, 1, just as one possible case. Then the 
flat directions tp and ip' are given by: 



with respective masses: 
(m 2 Ui + m\ + m\ 



(L{+L 2 2 +e 3 ) 
V3 



(12) 



(m 2 



Li 



(13) 

Note that <p and <p' together break all of the SM gauge 
symmetries, hence, they can affect thermalization of in- 
flaton decay products [J]. 

These two directions independently oscillate with am- 
plitudes ipo and ip' , respectively, and frequencies m v and 
m v i respectively. 

There are four light scalar degrees of freedom which 
correspond to tp and ip'. The total number of scalars is 
Ntotai = 28, with 18 contained in the u\d\di monomial 
and 10 in the L\Liei monomial. Since the SM gauge 
group is completely broken, we have ^broken — 12. Then 
Eq. |gj) results in N Hght = 28 - 24 = 4 . 

This implies that (p and ip 1 are the only light fields 
around. Their mass eigenstates and eigenvalues, given 
by Eqs. (fT2|) and (fl3)) . respectively, are constant in time 
and therefore not affected by their oscillations. There- 
fore, the condition in Eq. ([5]) will not be satisfied. There 
are time variations in the mass eigenstates and eigenval- 
ues of the heavy fields. However, as discussed in the pre- 
vious section, so long as cp and tp 1 have rotational motion, 
resonant particle production will be negligible [2(| Hl[ • 

There are many examples where the fields in the two 
monomials can have Yukawa couplings to each other. In 
these cases simultaneous F-flatness of the two directions 
is not guaranteed in general. For an example, let us 
consider the two directions represented by QQQL and 
LLddd monomials. Note that Q and d can be coupled to 
each other through quark Yukawas. If this is the case, 
then the two directions cannot develop large VEVs si- 
multaneously. 

However, Yukawas couple fields with the same color 
index to each other. One can then choose color indices 
such that Q and d from the same family have different 
colors. This implies that there is a subset of directions 
which are simultaneously D- and -F-flat, and hence, can 
develop large VEVs independently. The situation will be 
similar for the directions represented by udd and QLd 
monomials 7 . 



In the above mentioned cases, one finds Nu g ht > 4, 
which implies that there are light degrees of freedom be- 
sides the two oscillating flat directions. This however is 
not a sufficient condition to have preheating. To eluci- 
date this point, let us consider a toy model with a U(l) 
gauge group and four scalar fields (f>\, 4>[, 4>2, 4>2 with 
respective charges +1, +1,-1,-1. The corresponding D- 
term is 



V D ~ g 2 (l^il 2 - |0 2 | 2 + Wi\ 2 ~ \4>' 2 \ 
After the following field redefinition 



(14) 



X 



V2 

h - <h) 

V2 



X 



V2 
{<t>'i ~ 



V2 



we have 8 

V D ~ g 2 (<p'rx'r + f'ix'i + vrxr + fiXiY 



(15) 



(16) 



where R and I denote the real and imaginary components 
of a scalar field respectively. We consider flat directions 
with maximal rotation 9 : 

ipn = ^ocos(m p t + 6) , pi — ^osin(m v < + 8) 

ip' R = ip' Q cos(m v it) , (p'j = <p' BW.(m v it) , (17) 

where 9 is a constant phase. There is a total of 8 scalar 
degrees of freedom. Four of them tpn, ipi, <p' R , ip 1 , rep- 
resent the two complex flat directions. The mass matrix 
for the other four degrees of freedom can be diagonalized 
to obtain instantaneous mass eigenstates. It is readily 
seen from Eq. (|16[) that one degree of freedom 

Xi(t) = (<Po + ^o 2 ) _1/2 t PoCoa(™,< p it)^R + 'p' sm(m ip d)x'i 



+ (pocos(mtpt + 0)xr + (posin(m v t + 9)\i 



(18) 



has a mass ~ gi}p\ + ^o 2 ) 1 ^ 2 an d is superheavy. Spon- 
taneous breakdown of the U(l) by the VEV of two flat 
directions yields a Goldstone boson 

X2(t) = (^o + ^o 2 r 1/2 - ^o sin ("V%fl + ¥>o cos ("V*)X/ 



(,5osin(m v < + 9)xr + tpocos(m v t + 9)\i 



(19) 



7 The second version of Ref. [II also cites other examples of two 
independent flat directions such as LLe and QLd. However, si- 
multaneous .F-fiatness of these two directions requires that the 
L in QLd must belong to the same family as one of the two L's 
in LLe. Therefore one cannot have two independent directions 
parameterized as L = L = e = i/jq and Q = L = d = tp' (up to 
two independent phases). It is possible that Q, L, d, and e all 



have a non-zero VEV, but the situation will actually be similar 
to that of multiple flat directions, which will be discussed later. 
This is also the case for two directions represented by LLe and 
LLddd monomials. 

A general D-term has the same structure, see Eq. JS). It is 
therefore straightforward to generalize this discussion made for 
a U(l) gauge group. 

Generalization to the case with general rotation is straightfor- 
ward. 
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which is eaten up via the Higgs mechanism (hence re- 
moved from the spectrum). Note that in the special cases 
where cpo = or (p' — 0, the expression in Eq. (fT9|) is 
reduced to that for a single flat direction, given in the 
second line of Eq. ([7]) . 

Finally, there are two physical light degrees of freedom 



X3 W = (fo + Pa 2 ) 1 («lXfl + a 2X'i + a. 3 XR + a^Xi) 



and 
X*(*) = 

where 



and 



(Vo + ^o 2 ) 1 (biXR + hx'i + hXR + hX4 



a\ = —tpo 2 + ^o^o cos(Am< + 9) 

a-i — ~(po(p' sin(Ami + 6) 

a 3 = —(p' + ipofo cos(Am< + 9) 

o-4 — Vo^o sin(Ami + 9) , 

b\ = — ipQip' sin(Ami + 9), 

b 2 = (fia 2 - ipo(p' Q cos(Amt + 9) 

&3 = (fioVo sin(Ami + 9) 

bi = <p[) 2 — <po<p'o cos(Ami + 9) . 



(20) 



(21) 



(22) 



(23) 

These two eigenstates have a soft supersymmetry break- 
ing mass m x ~ O(TeV). Here Am = m v — m v ' is the 
mass difference of the two fiat directions. 

Time variation in Xs{t)i Xi(t) may lead to resonant 
particle production. Note however that this can hap- 
pen only if m v ^ m v '. Indeed it can be seen from 
Eqs. (I20l2ip that there will be no time-dependent terms 
m Xi(t) if the flat directions have the same mass. 

This can be qualitatively understood as following. The 
two rotating flat directions can be visualized as two vec- 
tors which rotate in the field space. So long as the flat 
directions have the same mass, the two vectors rotate 
with the same frequency, and so does the resultant vec- 
tor from their addition. Therefore the situation is effec- 
tively reduced to that of one rotating flat direction. In 
this explained in the previous section, all time 

variation is in the sector of heavy fields. Neither the mass 
eigenstates nor the mass eigenvalues of the physical light 
fields are time-dependent. In consequence, there will be 
no resonant particle production. 

Having Am 7^ is necessary but not sufficient for 
significant particle production. To demonstrate this let 
us consider two flat directions with hierarchically differ- 
ent VEVs. Without loss of generality, we assume that 
ip'o <C <po. Then, to the leading order in <p' /(po, we have 



Xs(t) 



-Xr 



(-) 



cos(Amt + 6)(x' R + XR) + sin(Ami + 6)(-x'i + Xi) 



and 

X4(t) = x'i 
sm(Amt + 9)(- X ' R + XR) 



(-) 

cos(Amt + 9)(x I + Xi) 

(25) 



The important point is that <p' /ipa appears in front of 
the time-dependent terms in Eqs. (|24I25|) 10 . Temporal 
derivatives of these terms which show up in the equa- 
tion of motion of X3(t), Xi{t) modes will therefore be 
suppressed by this factor. It follows from the analysis 
of [161 ] that modes with a physical momentum k < Ak 
will undergo resonance where 



(26) 



Had X3(t)> Xi{t) been massless (i.e. m x — 0), the width 
of the resonance band would have shrunk cx (ipo/ipo) 1 ^ 2 ■ 
Modes within the resonance band were exponentially am- 
plified, though the required time scale would increase by 
a logarithmic factor log((/?o / Vo) due to phase space sup- 
pression. However supersymmetry breaking inevitably 
results in a soft mass m x ~ O(TeV). Hence the reso- 
nance band completely disappears when 11 



(f' \ 77K 



(27) 



Since m x , m v , m^/ are all set by supersymmetry break- 
ing, even a small hierarchy between ip' and ipo will 
be sufficient to kill off the resonance band. Typically 
there will be no non-perturbative particle production if 
tpo > Wo 12 . 

We remind that the VEV of a flat direction at the onset 
of its oscillations is set by the higher-order superpoten- 
tial term of the form ip n /M n ~ 3 (induced by physics at a 
high scale M) which lifts it: </? ~ (m^M™" 3 ) 1 /™" 2 
Directions which are represented by different monomi- 
als usually are lifted by terms with different values of n, 
and hence have different VEVs. For example, QLd and 
QQQL are lifted by n = 4 terms, while udd and LLddd 
survives until n — 6 and n — 7 respectively [2] . Consider 
any two such flat directions ip' and ip which are lifted 
by superpotential terms of order m and n respectively 
(3 < to < n). For M = M P or M = A/ G ut, we indeed 



(24) 



10 In the opposite case when tpo C^J, the expressions for X3(*) an d 
Xi(t) are obtained by changing tp <-> ip' , Xr,i <-> /> (Ami + 
0) — » -(Amt + 8) in Eqs. (124125ft ■ The time-dependent terms 
are in this case suppressed by ipo/f'o <C 1. 

11 This point is ignored in the third version of [lTjl . see its added 
note. 

12 Note that Hubble expansion redshifts the VEV of two flat direc- 
tions the same. Hence Eq. H26|l determines the width of resonance 
band throughout the expansion of Universe. 
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have ipo > 10f' o for all cases with m ^ n. Then, as we 
explained above, there will be no preheating at all. 

One comment is in order before closing this subsec- 
tion. As we discussed in section 1, the masses induced 
by the VEV of flat directions are supersymmetry con- 
serving. This implies that the spectrum of fcrmions is 
the same as that of bosons except for the soft supersym- 
metry breaking effects. In particular, the eigenstates of 
fermionic partners of X3(t) an d Xi{t)-> denoted by X3(i) 
and Xi{t) respectively, have the same time dependence 
as m Eqs. (|24l25p 13 . Hence, similar to Eq. {26]), corre- 
sponding modes with a physical momentum k < Ak f erm 
have resonance instability where 



Ak 



ferm 



(— ) 



1/2, 



(28) 



The important difference is that fermionic modes have no 
supersymmetry breaking mass m x . Thus the resonance 
band shrinks as a result of a hierarchy ip' Q <C ipo but, 
unlike the bosonic case, does not completely disappear. 
Therefore one might think that flat directions will quickly 
decay via non-perturbative production of fermions. 

Let us estimate the rate for production of fermions 
The key point is that, due to Pauli blocking, the oc- 
cupation number of fermions cannot exceed 1. An oc- 
cupation number ~ 1 is achieved within a time 5t ~ 
Ak~jl rm [Til]. The energy density in the produced parti- 
cles is pferm ~ A{;J erm . Hubble expansion and rescatter- 
ings deplete the resonance band, after which production 
of fermions within the band will resume. Assuming that 
this happens immediately, an amount Afcj erm of energy 
density is transferred to fermions within each time in- 
terval ~ AkJ erm . Then the energy density in fermions 

increases by ~ (Afcy erm /if) in one Hubble time. Parti- 
cle production becomes efficient when this is comparable 
to the energy density in flat directions 14 . In a Universe 
dominated by relativistic particles (from inflaton decay) , 
the latter changes oc if 3 / 2 . With the help of Eq. (|2% 
it turns out that production of fermions will be efficient 
when if <C m v (m tp /ipo) 4 ^ 5 ■ This is sufficiently late in 
order for the flat directions to affect thermal history of 
the Universe O. 



B. Multiple directions represented by a polynomial 

When the family indices of lepton and quark multiplets 
are taken into account, we will have a gauge- invariant 
polynomial which represent multiple flat directions 15 . 



13 Note that we are only interested in the light fields. Hence we do 
not consider gauge and gaugino fields as they have a mass ~ gtpo 
and are superheavy. 

14 For hierarchical VEVs, i.e. ipo 2> p' , * ne ener gy density is 
mainly in tp, 

15 In this sense the H u Hd monomial is very special, since it is the 
only one which has no family indices. 



For example, let us consider the H u Li polynomial, where 
1 < i < 3 (for a detailed discussion, see [221]). 

Assuming D-flatness, one can always go to a basis 
where if* = L\ — L\ = L\ = (superscripts denote 
weak isospin components) , and (up to an overall phase) 



Hi 



Vo 
V2' 



rl _ ago rl _ bip rl _ \AH a l 2 -| b l 2 yo / 0Q \ 

L \ - -75 , L 2 - — ' ^3 -jn ' \ lJ ) 



V2 



\/2 



where |a| 2 + \b\ 2 < 1. 

A flat direction is then given by 



and its mass is: 



a*L\ 



b*L\ 



\b\ 2 H) 



V2 



(30) 



+ \a\ 2 m 2 Li + \b\ 2 m\ 2 + (1 - \a\ 2 - |6| 2 )r 



(31) 

Now the flat manifold has complex dimension 3. Note 
that Ntotai = 16 and Nb ro ken = 3, and hence Eq. ((4]) 
results in Nu g ht — 10. This implies that there are four 
light degrees of freedom besides those which parameterize 
the flat manifold. 

The trajectory of the flat direction VEV is determined 
by its initial value, the soft breaking terms and higher or- 
der superpotential terms which lift the flatness. The flat 
direction motion is confined to a complex plane when a, b 
are constant in time (in particular for a = b = 0, a = 
and \b\ — 1, or \a\ — 1 and 6 = 0). However, for a general 
motion on the flat manifold, a and b can assume any val- 
ues in the allowed ranges, and hence are time dependent. 
The pha se motion of the fields was already recognized in 
Ref. [H. 

However, note that just time variation of mass eigen- 
states is not a sufficient condition for a quick non- 
perturbative decay of the flat direction. A closer analysis 
suggests that the manifold spanning the flat direction is 
curved and the kinetic terms are non- minimal [22] ] . It was 
already noticed in Ref. [22| , that a small perturbation of 
the flat direction trajectory will be highly non-trivial, 
the trajectories have instability which leads to possibly a 
chaotic motion. Particle production due to a chaotic mo- 
tion is much more involved than that for an oscillatory 
motion and, a detailed investigation is needed to clarify 
this situation. 

Further note, that the other 4 light scalars (besides 
6, which parameterize the flat manifold) are expected to 
have (self-)interactions with gauge strengths (from the 
D-terms), as they do not correspond to genuinely flat 
directions, see Eq. (|29I30|) . Such interactions can st rong ly 
suppress non-perturbative particle production 0, HEf ■ 

An important point is to note that the flat manifold 
does not remain the same for all field values. As men- 
tioned, higher-order superpotential terms of the form, 
ip n /M n ~ 3 (with n > 3), lift the flatness for large field 
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values. If one writes all higher order terms which re- 
spect the SM gauge symmetries, then none of the MSSM 
flat directions survive beyond n — 9 2}. A higher or- 
der term ip n /M n ~ 3 becomes important at a field value 
\ip\ ~ (m ip M n - 3 ) 1 / n - 2 Therefore, terms with dif- 

ferent n take action at different VEVs and lifting of the 
flatness happens at stages. In some physically interesting 
cases the flat manifold has complex dimension "one" at 
very large field values, which effectively leaves us with a 
single flat direction. This happens for the LLe, LLddd 
and QuQue polynomials (when all families are taken into 
account). The flat manifolds corresponding to these poly- 
nomials have complex dimension 3,3, and 16 respectively. 
However beyond n — 4, equivalently at field values larger 
than (m^M) 1 / 2 , they all become effectively of complex 
dimension one [2| . These remaining single flat directions 
are eventually lifted with n = 6, n = 7, and n = 9 terms 
respectively [2j. Then, starting at large VEVs, the sit- 
uation will be similar to that for a single flat direction 
discussed in the previous section, and hence there will be 
no preheating. 

Finally, as a side comment, we emphasize that the 
flat directions can "terminate" an early stage of non- 
perturbative inflaton decay via preheating, regardless of 
their fate 

!• They can induce large masses to the in- 
flaton decay products at time scales long before the flat 
direction oscillations start. Hence, contrary to the stan- 
dard lore [HI , the inflaton can decay perturbatively even 
if it has large couplings to the matter fields. 

Within MSSM the only renormalizable interactions 
of a gauge singlet inflaton, <&, with matter can hap- 
pen through superpotential terms of the form h§H u L 
or h^H u Hd 0. Rest of the superpotential couplings 
to matter fields are through non-renormalizable terms. 
The couplings of a singlet to the matter fields through 
derivative terms (dictated by Kahler potential) will also 
be non-renormalizable. 

In models of large field inflation, a typical VEV of the 
inflaton is ($) ~ O(Mp), during and just after infla- 
tion. Then any renormalizable couplings of the infla- 



ton induces a large SUSY conserving mass, hMp, for H u 
and L, or H u and H e i multiplets. This implies that nei- 
ther H U L nor H u Hd directions can develop a large VEV 
during inflation. Hence neither of them are relevant in 
the process of thcrmalization, unless the inflaton is very 
weakly (i.e. gravitationally) coupled to all fields. 



IV. CONCLUSIONS 

To conclude, the particle production via preheating 
does not happen for a single flat direction, it is also not 
likely for multiple flat directions. In many cases which 
are of physical interest, after taking into account of the 
-F-term constraints and F'-term lifting, the multiple flat 
direction case becomes effectively the same as that of a 
single flat direction, or collection of independent single 
flat directions. In either cases preheating is not relevant. 
In particular, the flat direction which has the largest 
VEV survives long. The longevity of the flat direction 
then ensures that thcrmalization of the inflaton decay 
products is considerably modified 0, Q . 

The long lifetime of MSSM flat directions leads to late 
thcrmalization of inflaton decay products and a low re- 
heat temperature. This gives rise to a natural solution 
to the thermal overproduction of gravitinos, and other 
dangerous relics 0, Q • Regardless of their longevity, flat 
directions can also lead to a perturbative decay of a gauge 
singlet inflaton, even if it has large couplings to matter 
fields [![. To summarize, the MSSM flat directions in 
general play crucial role in reheating and thcrmalization. 
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